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Abstract 

We use a multi-scale similarity analysis which gives specific relations 
between the velocity, amplitude and width of localized solutions of non- 
linear differential equations, whose exact solutions are generally difficult 
to obtain. 

In this paper, wavelet-inspired approaches for localized solutions of NPDE 
are explored [1,2]. A first method provides relations between the characteristics 
of such solutions (amplitude, width and velocity) without the need of solving the 
corresponding NPDE. The method uses the multi-resolution analysis [2] instead 
of the traditional tools like the Fourier integrals or linear harmonic analysis 
which are inadequate for describing such systems. This scale approach has the 
advantage that it does not need the explicit form of the exact solutions. Hence, 
it is useful especially in situations when such solutions are unknown. The self- 
similar character of the fission process of fluid drops is an example where the 
same type of singularity occurs in any scale [3,4]. In the following we introduce 
a one-scale analysis (OSA) for the NPDE, in terms of their localized traveling 
solutions. The so called OSA analysis, described and applied here for localized 
traveling solutions belonging to any type of NPDE, provides algebraic connec- 
tions between the width L, amplitude A, and the velocity V of the solution, 
without actually solving the equation. The procedure consists in the substitu- 
tion of all the terms in the NPDE, according to the rules: 

i (substitution valid for traveling waves) Ut — > —Vux 

ii u—t +A, ~* ±A/L'^, 
and so forth for higher order of derivatives. This substitution in eq.(ii) is 

possible only for localized (finite extended support) solutions having at least 
one local maximum (like solitons or Gauss functions), if they exist. Since we 
are interested in traveling solutions, the first substitution reduces the number 



of variables from 2 to 1 so that wc arc now dealing with an ordinary differen- 
tial equation instead of a PDE. Then, the second substitution transforms the 
ordinary differential equation into an equation in the parameters describing the 
amplitude, width and velocity. Consequently, the NPDE is mapped into an 
algebraic equation in A, L and V. The proof of the method follows from the 
expansion of the soliton-like solution u{x, t) = u{s) withs = x — Vt, in a Gaus- 
sian family of wavelets ^'(s) = 7Ve'3(*), where Q{s) is a polynomial and the 
normalization constant [2]. If we choose Q = —is — {{s^2}\over2} we obtain a 
very particular wavelet with the support mainly confined in the (—1, 1) interval, 

namely 'J'(s) =exp —is — ^ j-n^l^. We have the discrete wavelet expansion of 
u 



k k 

in terms of integer translations fc of ^P, which provide the analysis of localiza- 
tion, and in terms of dyadic dilations 2-' of \1/ which provide the description 
of different scales. The idea of the proof is to choose a point where the ex- 
pansion in cq.(l) can be well approximated by one single scale such that 'u(s) 
can be approximated with a sum of phases in sthe OSA approach is given by 
the convective-dispersive equations, for which the most celebrated example is 
provided by the KdV equation 



Ut -I- uu^ -I- u^xx = 0, (2) 

In Table 1 we present examples of pure dispersive NPDE, identified in the first 
column by the form of the equation, and in the second column by a correspond- 
ing traveling localized solution, if the analytical form is available. Such exact 
solutions provide special relations between L, A and V , which are given in the 
third column of Table 1. In the last column we introduce the results of OSA, 
namely the relations between these three parameters, provided by eqs.(ii). The 
usefulness of the approach may be checked, by a quick comparison between the 
second and the third columns. The case of the KdV and MKdV equations, 
eq(2,3), are described in the first and second rows of the Table 1, and also are 
analyzed in previous papers of the same authors. Moreover, the same relations 
like in the KdV case remain valid even for the solutions of the "compacton" 
type [6] 



, , v^fccos Tfc (a; - 4A;^) 



(l - fcos k[x- AkHf 



where L = n/Gk, that is L~l/A, like in the Table 1. Next example (third row) 
is provided by a generalised KdV equation, in which the dispersion term is 
quadratic 

w* + {u^)x + (w^)xxx = 0. (4) 
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Eq.(3), known as K(2,2) equation because of the two quadratic terms, admits 
compact supported traveling solutions, named compactons [1,5,7-9]. In general, 
the compactons are obtained in the form of a power of some trigonometric 
function defined only on its half-period, and zero otherwise, in such a way that 
the solution is enough smooth for the NPDE in discussion. In the above example 
the square of the solution has to be continuous up to its third derivative with 
respect to x. Different from solitons, the compacton width is independent of the 
amplitude and this fact provides the special connection with the wavelet bases. 
The compactons are characterize by a unique scale, and it is this feature that 
makes it possible to introduce a nonlinear basis starting from a unique generic 
function. For eq.(3) the compacton solution is given by 



Vcix - Vt) = ^cos^ 



Vt 



(5) 



if Ix — yt| < 277 and otherwise. Here we notice that the velocity is proportional 
to the amplitude and the width of the wave is independent of the amplitude, 
L = 4. As a field of application wc mention that the quadratic dispersion term is 
characteristic for the dynamics of a chain with nonlinear coupling. The general 
compacton solution for eq.(3) is actually a "dilated" version of eq.(4). Actually, 
this combination is just a kink compacton joined smoothly with an antikink one 



Vkak{x-VT;X) = 



(6) 



AV 



■ cos 



x-Vt 



0..., 

for - 27r < x-Vt< 0, 



— for < x-Vt<X, 



■ cos 



Vt-X 



for A < X - Vt < X + 2iT, 
0... 



Finally, we can construct solutions by placing a compacton on the top of a 
KAK. Such a solution exists only for a short interval of time {X/V), since the 
two structures have different velocities. The analytic expression of the solution 
is given by 

yi{x,t) = 71kak{x - Vt; X) + [r,, {x - Vt - 27r) + — j ^ (7) 

for < t < (A — 47r)/(V^' — V) and zero in the rest. Here xi^) is the support 
function, equal with 1 for ja: < 1 and in the rest, and V' = 3max{77g/4} -|- 2V . 
For the K(2.2) compacton eq.(4) fulfills some relations between the parameters: 
A = AV/'S and L — 4 [7]. The relations provided by OSA in the last column 
of the third row, predict such relations, and hence also proove the existence of 
the compacton. Another good example of the predictive power of the method is 
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exemplified in the ease of a general convection-nonlinear dispersion equations, 
denoted by K(n,m) 



Vt + V:+vZ.^O, (8) 

Compacton solution for any n\neqm are not known in general, except for some 
particular cases. In this case we find a general relation among the parameters, 
for any n,m , shown in the fourth and fifth rows. These general relations L(A,V) 
approach the known relations for the exact solutions, in the particular cases like 
n=m (fourth row), n=m=2 (third row), n=m=3 (first reference in [1]) and n=3, 
m=2 ; n=2, m=3 (fifth row). These results can be used to predict the behavior 
of solutions for all values of n,m . 

The situation is different in the case of conipactons, which allow also sta- 
tionary solutions. When linear and nonlinear disspersion occur simultaneously, 
like in the so called K(2,l,2) equation 

Ut + ul+ Uxxx + e-ul^x = 0> (9) 

where e is a control parameter, the OSA yields a dependence of the form L = 
■\J^^^, which still provides a constant width if V = ±A + 2e . In this case, 
the speed is proportional to the amplitude, but can change its sign even at non- 
zero amplitude. Solutions with larger amplitude than a critical one {A^rit = 
=F2e)move to the right, solutions having the critical amplitude are at rest, and 
solutions smaller than the critical amplitude move to the left. This behavior was 
explored in [7], too. A compacton of amplitude A on the top of a infinite-length 
KAK solution of amplitude 5 

f X — Vt\ 

u{x,t) = Acos^ — - — + (5, (10) 



is still a solution of the K(2,2) equation, with the velocity given by ^ = 
I (2(5 + A). For A = —26 the solution becomes an anti-compacton moving 
together with the KAK. In the case of a slow-scale time-dependent amplitude 
the oscillations in amplitude can transform into oscillations in the velocity. The 
key to such a conversion of oscillations is the coupling between the traditional 
nonlinear picture (convection-dispersion-diffusion) and the typical Schrodinger 
terms. In Table 2 we present another class of NPDE, namely the dissipative 
ones. These equations generalize the linear wave equation (first row) where 
there is no typical length of the traveling solutions. The wavelet analysis pro- 
vides the correct expression for the dispersion relation {V = c ^ fc^ = uj'^/c^) 
with no constraint on either the amplitude A or on the width L . In the second 
row we introduce the Burgers equation which represents the simplest model for 
the convectivc-dissipative interaction. Dissipative systems are to a large extent 
indifferent to how they were initialized, and follow their own intrinsic dynamics. 
We provide in the second column an analytic solution of the Burgers equation. 
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For some special of values of the integration constants (2C < V^,D = )tlie 
solution becomes a traveling kink 

u{x, t) = V+ - 2C ta.n^^/V'^ - 2C{x - Vt)). (11) 

By applying the OSA approach to the Burger equation (third column) we obtain 
the same relation between amplitude and half- width, like in the case of the exact 
solution eq.(9), providing the velocity is proportional with the amplitude. In 
the following we apply the OSA approach to investigate a nonlinear Burgers 
equation 

Ut + auf-/xu^, + cuT = 0, (12) 

called quasi-linear parabolic equation [9] , and used to describe the flow of fluids 
in porous media or the transport of thermal cniergy in plasma. The existence and 
stability of waves or patterns is strongly dependent on the coefficients a, /i, c, /?, 
and 7 , and at this point the OSA can be useful again since there is no general 
analytic solution for eq.(lO). The result of the OSA approach is presented in the 
third row of Table 2. The typical scale of patterns depends on the parameters 
in the equations and the amplitude of the excitations, in a complicated way. 
However, in order to test OSA again, we found a simple class of exact solutions 
when c = , presented in the fourt row in Table 2, and expressed as the 
inverse; of a degenerated hypergeometric function. In this expression we have 
A= ijkV°'-'^F[{{k - l)a°')-\z] = ia/V)u"'-^ and a = {k - l)/(m - 1). The 
asymptotic behavior of the left hand side of the solution given in fourth row, 
second column, is described by 



r(a-l) 



(-1)°' -I 3_r(«-i)e« 



+ 0{l/z). (13) 



If z approaches + oo the solution increases indefinitely like an exponential. For 
a >1 (strong difussion efi'ects), for even k and for even (3 , the traveling wave 
u{x — Vt) has a negative singularity towards -oo at a; -|- xq = r(a + 1) (— 1)" < 
. For k odd there is also a singularity at a; -|- a;o > . If A: is even and /3 is 
odd (the singularity is pushed towards imaginary a; ), or if < a < 1 , the 
singularity is eliminated and the solution becomes semi-bounded, like in the 
particular situations investigated in the article [9]. In this case, OSA provides 
again the correct relations, since we obtaine the special behavior of the solution 
if the velocity is proportional to the power /? — 1 of the amplitude A . Also, 
we predict the space scale of these semi-compact pulses, namely the length 
^ ~ v±a/3A'^-^ • OSA analysis can be applied in the case of sine-Gordon 
equation, fifth row of Table 2. The solutions with the velocity proportional with 
are characterized through the OSA approach by a transcendental equation 
in A, identical with the equation fulfilled by the amplitude A of the exact sine- 
Gordon soliton. In the sixth row, we present the cubic nonlinear Schrodinger 
equation (NLS3) which has a soliton solution. In the sixth row of Table 2 we 
present the NLS3 equation together with its one soliton solution of amplitude tjq 
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, obtained by the inverse scattering method. In the last cohimn we also show the 
relation between the parameters of a localized solution, obtained by OSA. The 
equation for L{A, V) is more general than that one fulfilled by the soliton, and 
hence is related to more general localized solutions. By choosing the velocity 
proportional to the amplitude, we reobtain the L ^ 1/A = I/tiq ^ 1/V typical 
relations for the soliton given in the second column. 

"2^^*"" + ^^"^^* + ''*'"°' ^^^^ 

then the L parameter gives an estimation for the wavelength of the wavefunc- 
tions, or for the correlation length in a Bose model 

L = ^ ~ ^ — for small A. (15) 

^y2m{E -V) + aA^ ^1m{E - V) 

For the general case of a NLS equation of order n (seventh row), where a general 
analytical solution is unknown, the method predicts a special L = I/(^, V) 
dependence, shown in the third column. Contrary to third order NLS, where the 
dependence of L with A is monotonous for V ±A (n = 3), at higher orders 
than 3, the L{A) function has discontinuities in the first derivative. This wiggle 
of the function n = 4 hold at a critical width, possibly producing bifurcations in 
the solutions and scales. As a consequence, initial data close to this width can 
split into doublet (or even triplet, for higher order NLS) solutions, with different 
amplitudes. Such phenomena have been put into evidence in several numerical 
experiments for quintic nonlinear equations [8-10]. The final example of Table 
2 is provided by the Gross-Pitaevski (GP) mean field equation, which is used to 
describe the dilute Bose condensate [11]. The scalar field (or order parameter) 
governed by this equation was shown to behave in a particle manner, too, since 
it can contain topological deffects, namely dark solitons. The space scale L of 
such solutions is important, for both the theory and experiment, since is related 
to the trap dimensions and to the scattering length. In the last row of the Table 
2 we give one particular solution of a simplified one-dimensional version of the 
GP equation [12] 

\ ' = A + v;xt(^) + ^;^*(^' t) ) *(^' (16) 

where a is the s-wave scattering length and Vg^t is the confining potential. In the 
solution provided in the table, the half- width of the exact nonstationary solution 
is L = 1/ \/i>\— p^, where Vc = \/l — aV is the Landau critical velocity, and p = 
dq{t)/dt is the momentum associated with the motion of this disturbance. It is 
easy to check that the OSA provides a good match with this exact solution, and 
also L fits the correlation length Iq = ^ ^J^-i^ ■ We stress that such estimation 
of the length is also important in nuclear physics where one can explain the 
fragmentation process as a bosonization in a-particles, inside the nucleus. Such 
systems are coherent if the wavelength associated with the cluster (the resulting 
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L in the GP equation) is comparable with the distance between the a -clusters. 

ut + au^ + nul^ + c<^^ = 0. (17) 

Here to, k and n are integers and the corresponding terms are responsible of 
the nonlinear interaction (convective term), dissipation and dispersion [9]. The 
above equation is related to weakly nonlinear phenomena, and it occurs in mod- 
eling porous medium, magma, interfacial phenomena in fluids (and hence ap- 
plications to drop physics), etc. Thi OS A approach maps this equation into 

aTOA™-^ - VL^ - fik'^A'^-^L + n^A^'^ = 0, (18) 

Table 3, first row. The most symmetric case is obtained when either V = 

(stationary patterns) oi V ^ A™^^ . In this situation the condition to have a 
monotonous dependence of L as a function of ^ is 2fc = m -|- n which yields a 
scale structure 



[ 2m{a-Vo) )^ ' ^^^^ 

where we put V = mVoA"^~^ . The condition 2fc = to + n is just the condition 
obtained in [9] from a scaling approach. This condition assures the universality 
of the corresponding patterns, and it is the unique case in which L depends on 
a power of A. In the above cited paper, the author finds out the condition for 
mass invariance at scaling transformations asTO = n-|-2 = A;-|-l. In our case 
we just have to request the product AL (which gives a measure of the mass, 
or volume of the pattern, like in the case of one-dimensional solitons) to be a 
constant. This gives the condition fc — to -|- 1 = which, together with the 
general invariance condition 2k = m + n , reproduces TO = n-|-2 = A;-|-l. In 
this case we have patterns characterized by a width 



_ ± yZ/i^fc'* - 4TO,n3(a - Vb) 

2mA{a-Vo) ^A^ik^' ^ ' 

If a ~ Vo the width approaches n^/A^k"^ . In order to make L independen of 
A , like in the compacton case, we need m = k, which together with the first 
invariance condition 2k = m + n, yields m = n = k. This is the exceptional case 
when the dissipative and dispersive processes have the same scaling, resulting 
form the invariance of the eq.(12) under the group of scales. Finaly, if we choose 
L ^ V we obtain the condition A; + 1 = 2m which (together with 2k = m + n) 
is the condition for spiral symmetry and occurence of similarity structures [9]. 
The next example is provided by one of the most generalized KdV equation, 
which is generated from the Lagrangian [5] 



C{n,l,m,p) = j 



2 (p + l)(p + 2) — 2 



dx, 
(21) 
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where a, f3 and 7 are parameters adjusting the relative strcntgh of the inter- 
actions, and n, I, m, p are integers. For example, for 7 = one re-obtains 
the K(2,2) equation, and for 7 = m = 0,p = 1 one obtains the KdV equa- 
tion. The associated Euler-Lagrange equation in the function (p^ = u{x, t) — > 
u{x — Vt) = u{y), reads after one integration 

= ^""^^ - Pmu^-'UI + 2I3{VJ-Uy)y + ^-^U^-\Uy)\Uyyf (22) 
7Z 

- 2-(«"K)'"^Ky)')y + 7(w"(%)'uyy)TO + C, 

where C is the integration constant. By using the OSA we obtain the following 
important result, expressed in the second row of Table 3: The unique case 
when such an equation allows compact supported traveling solutions is when 
m = p = n + r,C = Q and V = VqA"'- . This result is in full agrement 
with the variational calculation in [5]. Both eqs.(12) and (16) are rather more 
qualitative than capable of modeling measurable phenomena. That is why we 
introduce now a more general model equation, in the form 

Ut + f{u)x + g{u)xx + h{u)xxx = 0, (23) 

where /, g and h are differentiable functions of the the function u{x,t) itself. 
The OSA approach gives the equation 

Ag"{A)+g'{A) A^ h'" {A) + 3Ah" (A) + h' (A) 

A general analysis of eq.(17) is difficult, and the best ways are numerical in- 
vestigations obtained for particular choices of the three functions. Wc confine 
ourselves here only to show that the class of solutions which have similarity 
properties are those for which V = Vof'{A) . In this case eq.(18) can be re- 
duced to 

L2/'(1 - Vo) + L{Ag" + g') + A^h'" + ^Ah" + /i' = 0, (25) 

case which is presented in the third row of Table 3. This last relation can be 
used for different purposes. For example, given a certain type of dispersion 
and difusion { g,h fixed), we can estimate for what types of nonlinearity ( /) 
the width L will have a given dependence with A . Or, if we know for instance 
f{u) = /ou'^i and h{u) = Hqu'^^ , we can ask what type of diffusion g we need, to 
have constant scale (width) of the patterns (waves), no matter of the magnitude 
of the amplitude A . In other words, which is the compatible diffusion term, for 
given nonlinearity-dispersion terms, which provides fixed scale solutions. The 
result is obtained by integration eq.(19) with respect to g{u) 

gin) = (l + g2 + u^^ - ^:ML_^„9. + CsLog u + C, (26) 

where C-i,i are constants of integration. In a similar way one can check the 
existence of different other configurations by solving eq.(19), or more general. 



-V + f'iA) 



8 



eq.(18). A last application of this method, occurs if the KdV equation has an 
additional term depending on the square of the curvature 



+ Wx + Uxxx + e (ulx)^ = 0. 



(27) 



This is the case for extremely sharp surfaces (surface waves in solids or granu- 
lar materials) when the hydrodynamic surface pressure cannot be linearized in 
curvature. Such a new term yields a new type of localized solution fulfilling the 
relations 



If we look for a constant half- width solution (compacton of 1/L = a) we need a 
dependence of velocity of the form V {l + a'^e/8)A + l/8eA + a/4 . There are 
many new effects in this situation. The non-monoton dependence of the speed 
on A introduces again bifurcations of a unique pulse in dublets and triplets. 
Also, there is an upper bound for the amplitude at some critical values of the 
width. Pulses narrower than this critical width drop to zero. Such bumps can 
exist in pairs of identical amplitude at different widths. They may be related 
with the recent observed " oscillations" in granular materials [7] . The examples 
presented in Tables 1-3 prove that the above method provides a reliable criterion 
for finding compact suported solutions. 
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Table 1: Traveling localized solutions for nonlinear dispersive equations. 



NPDE 


Analytic solution and the 




OSA 




relations among parameters 




approach 


Ut + 6uUx + Uxxx = 


Asech^^-/*; L = 


L = 


|y ±6A|-V2 




V = 2A 


IfVr 




Ut + V^Ux + Uxxx = 


ylsech^-/*; L = l/A, 


L = 


|y±6A2|-l/2 




A = VV 


If y 





/ \ 1/2 

Ut + iu% + {u^)xxx = Acos^^, if \{x-Vt)/4\<n/2; L = / ^||^ j 

L=4 



Acos 



2 1 x-Vt 



if la; - Vt\ < 



2n7r 
n-1 



1/2 



and else; 

r 4n A 2Vn 



if y = 



Ut + + = 

n ^ m 



unknown 
in general 



^ _ n(n" + l)A"~^ 



1/2 
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Table 2: Traveling localized solutions for nonlinear diffusive equations. 



NPDE 


Analytic solution and the 
relations among parameters 


OSA 
approach 


Utt — \-utt = 


fc2 = wVc^ 


V = c 
A, L arbitrary 


Ut + UUx - Uxx = 


VC y2tan(VC 1/2:1.-/* +£») 

+v 


L = {A±V)-'^ 
If y ~ A, L~ 1/A 


Ut + a{u'")x - f^{u'')xx + cu'^ 
=0 


only particular cases 
known 


cA-^L^ + {V±amA"'-^)L 
±fik'^A''-'^ = 


Ut + o(u™)x - fJ,{u'')xx = 


-Az"" iFi{a,a+ 1, z) = x + xq 


am — a ' 

if = 


Uxt — sin u = 


1 x-Vt 

A tan 7 6 J- 


±^^ = sin A 
IfV = L^,A = sinA 


+ + 21^-12* = 


ryoe'('^*+'=^)sec/i[r?o(a; - Vt)]; 
L = 1/% 


J ±v±^|y2-4A2| 
IfA^V, L = l/A 


+ *^^ + |^|""'* = 


unknown in general 




i^t = 
+ [a|*|2 + - 1]* 


ip+ \Jv1^ — p2 X 
tanli[a-\/?;2 — p2(^2;_ 


L = (aA2±T/-l)-i/2 

If y ~ ±i,L~ i/(Av^) 
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Table 3: Traveling localized solutions for dissipative-dispersive equations. 



The NPDE equation 



OSA approach 



T _ Am-k fJ.k'^±\/t^'^k^-4:mn'^{a-Vo) 
^ — ^ • 2m(a-Vo) 

if V = mVoA''-'^ 



+ l^u'-HUyY{Uyyf - {Uyr^Uyyf) y 

+ 7(M"(7i,,)'M,,y),,y +C 



2L'+4((n + / + l)Vo - a) 

-2L'+2(/ + n + l)(/ + n + 2)/3 
-(/ + n + 1)(2 + 2n2 + 3/ + f + + :il))-f 
if C = 0. \' = lo-l'" iiiid = p = II + I 



Ut + f{u)x + 9{u)xx + h{u)xxx = 



g' + Ag"T[{.Ag" + g'f 

-4/'(l - Vo){A^h"' + 2Ah" + h') 

x{2A'^h"' + 6Ah" + 2h')-'^ 
iiV = Vof'iA) 



111- 
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